Potential Energy Landscape Equation of State 
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Depth, number, and shape of the basins of the potential energy landscape are the key ingredi- 
ents of the inherent structure thermodynamic formalism introduced by Stillinger and Weber [F. H. 
Stillinger and T. A. Weber, Phys. Rev. A 25, 978 (1982)]. Within this formalism, an equation of 
state based only on the volume dependence of these landscape properties is derived. Vibrational 
and configurational contributions to pressure are sorted out in a transparent way. Predictions are 
successfully compared with data from extensive molecular dynamics simulations of a simple model 
for the fragile liquid orthoterphenyl. 

PACS numbers: 64.70.Pf, 61.20.Ja, 61.20.Lc 



Recent years have seen a resurgence in studies de- 
voted to modeling the thermodynamics of supercooled 
liquids p], p, 0, 0, 9. Such studies aim to elucidate the 
physics of the liquid-glass transition, to develop a ther- 
modynamic description of out of equilibrium systems and 
to provide keys for a deeper understanding of the dynam- 
ics of supercooled states ||. Numerical studies, boosted 
by increased computational resources which now allow 
simulations to track the slowing down of the dynamics 
over more than 6 decades in time, are providing quantita- 
tive estimates for the free energy of simple model systems 

H pi [To|. The availability of such data provides strin- 
gent tests of the theoretical predictions |j| [l0|, |ll| and 
helps in the understanding of basic mechanisms associ- 
ated with the behavior of thermodynamic and dynamic 
quantities close to the glass transition. 

Among the thermodynamic formalisms amenable to 
numerical investigation, a central role is played by the 
Inherent Structure (IS) formalism introduced by Still- 
inger and Weber Jl^. Properties of the potential en- 
ergy landscape (PEL), such as depth, number and shape 
of the basins of the potential energy surface are calcu- 
lated and used in the evaluation of the liquid free en- 
ergy H 10. 13 1 In the IS formalism, the system free 
energy is expressed as a sum of an entropic contribution, 
accounting for the number of the available basins, and 
a vibrational contribution, expressing the free energy of 
the system when constrained in one of the basins [ fl2| . 

Important progress has been made after the discov- 
ery that — for models of fragile liquids — the number 
^(e/s) of distinct basins of depth ejs in a system of N 
atoms or molecules is well described by a Gaussian dis- 
tribution 0, 
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have also shown that the basin free energy can be written 
as the depth eis plus a vibrational contribution which, 
in the harmonic approximation, has the well known form 

M 

F mb (e IS ,T) = fcaT^lnGSfe^e/s)), (2) 

i=l 

where (jj(ers) is the i-th normal mode frequency (i — 
1...M) and j3 = l/k B T. The M normal mode frequen- 
cies define the shape of the basin. If relevant, anharmonic 
corrections can also be accounted for jlll [l3| . The quan- 
tity X)j=i ^ n i UJ i{ e is)/ UJ o) (where uj is the frequency unit) 
is found to depend linearly on the basin depth |L0| , i.e. 
can be written, in terms of two parameters a and 6, as 



M 



\]jx(u>i(e.is)/u ) = a + beis. (3) 

Hence, the vibrational free energy can be written as 

F mb (e IS ,T) = F mb (E ,T) + k B Tb(e IS - E ). (4) 

Within the two assumptions of Eq. (|l|) — Gaussian dis- 
tribution of basin depths — and Eq. (||) — linear de- 
pendence of the basin free energy on e/s — an exact 
evaluation of the partition function can be carried out. 
The corresponding Helmholtz free energy is given by |lC| ] 



F(T) = -TS conf (T) + (e IS (T)) 

+ F mb (E ,T) + k B Tb((e IS (T)) 



E n 



where 



(eis(T)) = (E - ba 2 ) - /3a 2 = eoo - /3a 2 



and 



Here the amplitude e accounts for the total number 
of basins. Numerical studies of models for fragile liquids 



S CO nf(T)/k B = aN -((e IS (T))-E ) 2 /2a 2 
= S x /k B -bpcr 2 ~(3 2 a 2 /2. 



(5) 
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FIG. 1: T dependence of (a) S con f, (b) ~£i H"i/"o)/N, 
and (c) (e/s) (per molecule) at the five studied volumes Vk 
(symbols) . Data are from Ref. jl3| . The full curves are simul- 
taneous fits of the three sets of data according to Eqs. (m), 
(^) , and ([]) . The dashed curves are the constant P paths (at 
P = and 200 MPa) calculated according to the IS-EOS as 
discussed in the text. The frequency unit is u}q = lcm~ 
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FIG. 2: V— dependence of the three combinations of PEL 
parameters contributing to the linear (squares), constant (cir- 
cles) and T -1 (triangles) components of P, according to 
Eq. ([ll]). The solid lines are the polynomial fit used to eval- 
uate Pconst, P T and Pi/t- Tables of the coefficients of the 
polynomial fit will be reported in Ref. Ju|. Here eoo is ex- 
pressed in 10 8 J /mol, Sac — fcs(a + be W 6 J/mol/K, and 
a 2 /2k B in 10 6 JK/mol. 



determined by the V dependence of the landscape prop- 
erties a, a, E , a and b. Eq. (|) shows that P can be 
split into three main contributions: a configurational one, 
Pmnf — related to the change in the number of available 
basins with V; an eis one, P eis — related to the change 
in basin depth with V; and a vibrational one, P v n, - 
related to the change in the shape of the explored basin 
with V. The T dependence of each contribution can then 
be studied independently. The explicit expressions for 
these contributions are: 
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con} (T, V) — T — Soc ( 



In the above equations, eoo and Soo are defined as the 
value of (ejs) and S con f at infinite T. Eqs. (^-(Q) show 
that, along constant volume V paths, the behavior of the 
thermodynamic quantities is controlled by the values of 
the PEL properties, as given by a, a, E (from Eq. (Q)) 
and by a and b (from Eq. (0)). 

In this Letter we study the volume dependence of 
Eq. (||) to provide an expression for the equation of state 
(EOS) based completely on landscape properties (PEL- 
EOS) This study provides a significant insight into 
the understanding of the pressure P and opens the way 
for detailed comparisons between experimental measure- 
ments (usually performed at constant P) and theoretical 
approaches based on the IS formalism. It may also help 
in developing an IS-based thermodynamic description of 
out-of-cquilibrium (glass) states and a theoretical defini- 
tion of the concepts of Active P and Tp6[. 

In thermodynamics, P is defined as the (negative) V 
derivative of the Helmholtz free energy. Hence P if fully 



P eis (T,V) = -A eoo + ^(a 2 /k B ) (9) 



Pvib(T, V) = -k B T—{a + b eoo ) + g^(bo- 2 ). (10) 

By grouping together all the contributions with the 
same T dependence, P can be expressed in term of V 
derivatives of only three combinations of the five PEL 
parameters O] 



where P, 



p{t, v) = p const + tp^ + r- x p 

i/T 0V y 



(11) 



-de^/OV, P = dSoo/dV - k B 



& eoo ) and P i/t = ^(a 2 /2fc B ). 

The present formalism also provides an expression for 
the so-called inherent structure equation of state (IS- 
EOS), Pis(Tis,V) H H, [20[ ||, i.e. the relation be- 
tween the pressure and volume of the typical IS and the 
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FIG. 3: Comparison between P evaluated according to the 
V derivative of PEL properties (lines) and P evaluated in the 
MD simulation using the virial expression (symbols). Solid 
symbols are equilibrium values. For each Vk, the symbol * 
indicates Pis for the coldest equilibrated state point. Open 
symbols are MD data calculated during a constant heating 
procedure starting from the IS configuration marked with *, 
as explained in the text. Lines are PEL-EOS for equilibrium 
(full lines, Eq. (pd])) and for the heating procedure (dashed 
lines, Eq. @). 



temperature T/<j of the equilibrium ensemble from which 
configurations were extracted. Indeed, the constant V 
steepest descent minimization procedure, which realizes 
the search for the closest local minimum starting from an 
equilibrium configuration, suppresses all the vibrational 
components (hence P v a> = 0) but keeps P CO nf and P eis 
frozen at their "equilibrium" initial value. As a result 
Pis, a purely mechanical quantity, can be expressed as 



-Bv Eo + Tls W So ° 



Pis{Tis,V) = P conf + 
1 d 



.2 



k B Tis dV 



K/2). (12) 



The present approach also predicts the behavior of P 
when an IS configuration is heated from T — at con- 
stant V. While the system remains in the same basin, P 
is given by Pis(T IS , V) + P vib (T, T I3 , V) where the only 
(linear) T-dependent part arises from 
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FIG. 4: Total (P), vibrational (P V ib) and eis plus configu- 
rational (P eis + Peon/) contributions to P along constant V 
(dashed lines) and constant P (solid lines) paths. At constant 
V, P vib is linear in T (Eq. @). 



We now apply the present derivation to the simple 
Lewis and Wanstrom (LW) model for the fragile molecu- 
lar liquid orthoterphenyl (oTP) (l3j |§. The LW model 
is a rigid three-site model, 



Pvib(T,Tis,V) = -k B T—(a + beis(Ti S )). 
ov 



(13) 



with intcrmolecular site-site 
interactions described by the Lennard Jones (LJ) poten- 
tial p2| . Its simplicity allows one to reach simulation 
times of the order of fis |l3| . In this model, as in the L J 
case, the anharmonic contributions are negligible, ejs(T) 

goes as 1/T, and J2i=i ^ n ( ljJ i( e is)/ UJ o) is linear in erg jl3). 
Hence this model is a perfect candidate for testing the 
validity of the PEL-EOS proposed here. 

We use the excellent data base of state points presented 
in Ref. |p| i) to calculate the V dependence of the PEL 
parameters; ii) to derive the EOS for the oTP model, 
and Hi) to compare it with the "exact" EOS calculated 
using the virial expression, as commonly implemented in 
molecular dynamics (MD) codes. 

Fig. |l| shows, for five constant V paths, the simulta- 
neous fit of the T dependence of S con f, the basin cur- 
vatures, and of (eis{T)), according to Eqs. (0), @, and 
(^). The possibility of fitting simultaneously, with the 
same values of a, a, Eq, a and b, the quantities (e/s(T)), 
J2i=i l n (^i(e/s)/w ) and S con f(T), supports the validity 
of the two main assumptions, i.e. Eq. (|l|) and Eq. (|4|). 

The V dependence of the three combinations of fit pa- 
rameters, — eoo, — fcs(a + be 00 ) and a 2 /2ks is shown 
m Fig.|. P{V,T) can be calculated from the V derivative 
of these quantities according to Eq. dill) and compared 
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FIG. 5: Inherent structure equation of state. Symbols are 
MD calculations, lines are PEL predictions. 



with the MD data. Such comparison is reported in Figj^. 
The very good agreement between the two set of data at 
all V's confirms that a very accurate EOS based on PES 
properties has been derived for this oTP model. 

The T dependence of the individual contributions can 
be evaluated according to Eqs. (|8|)-(^0|). Fig. || shows 
- along a constant V and a constant P path — P„j& 
and P con f + P eis - We note that at constant V, both P 
components are increasing with T, while at constant P, 
Pconf + Pe IS decreases on heating to compensate for the 
increase of P v a- 



The PES-EOS allows us also to contrast the iso- 
baric and isochoric T dependence of S con f, (e/s(T)) and 
12iLi m (wi/w )- Such comparison is reported in Fig. [l]. 
Here we note the faster decrease of S con f(T) and (e/s(T)) 
along constant P paths as well as the different trend in 
the change of basin shape |f23f. 

As a further check of the quality of the calculated EOS 
for the oTP model, Fig. [| compares the MD (virial) and 
PEL (Eq. (U)) IS-EOS. Since for the oTP model the 
term linear in T in Eq. ( |l2] ) is negligible |T?|] , Pjs fol- 
lows, to a good extent, a 1/Tjs law. The quality of the 
comparison supports the interpretation of P/g as the V 
derivative of the depth and number of basins sampled in 
the corresponding thermodynamic equilibrium state. 

Finally, Fig. | compares MD data (open symbols) and 
PEL (dashed lines) predictions in a run where T is in- 
creased starting from the T = IS configuration, as pre- 
viously discussed. The entire simulation is much shorter 
than the time needed to change basin, so that only the 
vibrational degrees of freedom are thermalized. Also in 
this case, the PEL expression accounts for the observed 
linear increase of P. 

The present EOS, based exclusively on PEL proper- 
ties, can and will be used to address important issues 
in the thermodynamics of supercooled liquids ^4j, such 
as the study of the intrinsic limit of stability of the liq- 
uid state |l5|, [2(], [HI a nd the IS-based thermodynamic 
description of aging|16| . 
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